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We develop exact, simple closed form expressions for partition functions associated with relativis- 
tic bosons and fermions in odd spatial dimensions. These expressions, valid at high temperature, 
include the effects of a non-trivial Polyakov loop and generalize well-known high temperature ex- 
pansions. The key technical point is the proof of a set of Bessel function identities which resum low 
temperature expansions into high temperature expansions. The complete expressions for these parti- 
tion functions can be used to obtain one-loop finite temperature contributions to effective potentials, 
and thus free energies and pressures. 
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O ■ I. INTRODUCTION 
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The applications of finite temperature field theory are numerous and diverse For many applications, a high- 
temperature expansion of one-loop contributions to thermodynamic functions is necessary. A typical one-loop con- 
tribution to a (d + l)-dimensional effective potential from a bosonic degree of freedom and its antiparticle has the 
] form 
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where the relativistic energy ujk is given by y^ + JP. In many cases, the mass M is a function of other quantities, 
most notably vacuum expectation values of fields. When the total effective potential attains its minimum, V may be 
identified as a contribution to the total free energy, and —V as a contribution to the pressure. A more general case 
is obtained when there is a non-trivial, but spatially uniform, Polyakov loop as well as a non-zero chemical potential 
/x.In this case we have 
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Note that the effect of a non-trivial Polyakov loop is to add a phase factor exp (±i9) to exp (— /3a;). We will evaluate 
Vb for arbitrary 9. In principle, the effect of /j, can be included by a careful analytic continuation 9^9 — but 
we do not consider it here. 

We will develop a high temperature expansion for Vb (9), valid for dodd, as well as a similar expression for the 
corresponding fermionic quantity Vp(9). This derivation is simple and exact, and generalizes the results of Dolan 
and Jackiw (J], who gave approximate high-temperature expressions for V B (9 — 0) and Vf{9 = 0) valid up to order 
M i ln(PM) in four dimensions. The work of Dolan and Jackiw was extended by Haber and Weldon || [|| who gave a 
complete expression for the Bosonic case Vb as an infinite sum over hypergeometric functions. Their work included a 
non-zero chemical potential. Later work by Actor showed that similar high-temperature expansions could be obtained 
using zeta-function techniques QS]. In both cases, higher-order correction terms are given by infinite series in f3M. 
Our expressions effectively resum these corrections into a simpler form. Analytical results for the case of a non-trivial 
Polyakov loop, 9^0, were first given in the case M = by Gross, Pisarski and Yaffe and by Weiss flT| . Our 
work generalizes their results to the case M ^ 0. The higher-order terms in our expressions are manifestly periodic 
in 9. This periodicity is important in the application of these results to our recent work with Miller on models of the 
deconfinement transition Jl^ . In this work, the eigenvalues of the Polyakov loop serve as the order parameters for 
deconfinement, a point of view also emphasized recently by Pisarski |13| . 

Before beginning the derivation, we give some examples of its application. As a first example, consider a scalar 
boson in the fundamental representation of an SU(N) gauge group. The Polyakov loop is an N x N unitary matrix 
given in general by 
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where T on the right-hand side indicates Euclidean time ordering. Here we assume that the Polyakov loop can be 
made spatially uniform by an appropriate choice of gauge. A global unitary transformation then puts P into the 
diagonal form 

Pjk = Sj k cxp (iOj) (4) 

and the partition function is 

X)Vb(fi). (5) 
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As a second example, consider the case of the gauge bosons themselves, which lie in the adjoint representation of 
the gauge group. The Polyakov loop in the adjoint representation is an (N 2 — 1J X (iV 2 — l) matrix. The partition 
function for the N 2 — 1 particles is 
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where the Sjk removes a singlet contribution, and the factor of 1/2 corrects for overcounting since Vb has both a 
particle and antiparticle contribution. The factor s accounts for spin degeneracy; in 3 + 1 dimensions s — 2, a 
consequence of the two possible polarization states of gauge bosons . 

For our third and final example, consider the evaluation of fermionic partition functions, which can be reduced to 
the general bosonic problem. A typical fermionic contribution of particle and antiparticle has the form 

V F (9 - z/3/i) = -\ [ In [1 + e-^+M+m _ I / m [l + e -^ h -ie-M (7) 
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which is easily written as 

V F (6) = -V B (it + 9). (8) 
For fermions in the fundamental representation of SU (N) , the partition function is 

«£>F(0i) = -«£M*-+0i) (9) 
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where the factor s again accounts for spin degeneracy. 

In section 2, we review the derivation of low temperature expansions for Vb (0) and Vp (0)- Section 3 derives the 
Bessel function identities which convert these low temperature expansions to high temperature expansions. Section 4 
applies these identities to the case of three spatial dimensions. A final section gives brief conclusions. There are two 
appendices. 



II. LOW TEMPERATURE EXPANSION IN D DIMENSIONS 

A low-temperature expansion for Vb (0) can be generated for arbitrary spatial dimension d by expanding the 
logarithm and integrating term by term, first over the surface of a c£-dimensional sphere, and then over a radial degree 
of freedom k [Mpl. The result, given in terms of modified Bessel functions [^4|, is 
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which is derived in detail in Appendix A. Each term in the series represents the contribution of n particles or antipar- 
ticles, with a corresponding phase factor of exp(±zn#). If the one-loop finite temperature functional determinant is 
represented as a functional integral over a space-time variable x^, the phase factors are associated with paths which 
wind non-trivially in the Euclidean time direction. 
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For fermions, we have 
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In a path integral representation, the factors of (—1)™ are a consequence of fermionic antiperiodic boundary conditions. 
We next derive a set of identities which resum these series for d odd. 



III. BESSEL FUNCTION IDENTITIES 



We will derive a set of identities for sums of the form 
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for p an even integer. Our starting point is the identity 
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which may be found in |jl4|| ; we provide a derivation in Appendix B which provides some physical insight into its 

origin. The notation is used to indicate that singular terms, here the 1/ |/| term, are omitted when I = 0. 
Using the recursion formula 



it follows immediately that 
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This in turn implies that 
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where C(^>)is an unknown function to be determined. Integration yields immediately 
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The function C(0) is determined from the behavior of K v {z) for z — > 
in combination with the standard result |l4fl 
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which is valid for < <j> < 2n. The right hand side of Eq. (|20j) is a rescaling of the second Bernoulli polynomial; it 
can be extended to all real values if </>is replaced by <f>mod2iT on the right hand side of the equation. This implies the 
leading behavior of the sum as z — > is given by 
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giving us finally 
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where we have introduced the notation </> + to represent fymodlit. When discussing fcrmions, we will also use 0_ to 
similarly represent an angle chosen to lie between — ir and tt. Note that the last part of this expression is automatically 
periodic due to the sum over I. 

Application of this technique a second time gives 
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which is needed for the case d = 3. We have used the standard result DA 
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Formulas appropriate for d = 5, 7, .. can also be derived in the same manner. 
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IV. HIGH-TEMPERATURE EXPANSIONS FOR D = 3 



We now can write complete expressions for Vb {&) and Vf (6) in three spatial dimensions: 
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Parts of this complete expression have been known for some time. For 8 = 0, the leading behavior is 



V B (9 = 0) 



M 1 



90/3 4 12/3 2 



M 3 
fa/3 



M 4 

167^ 



7 "4 



(25) 



(26) 



an expression first derived by Dolan and Jackiw [Q. The leading T 4 behavior for 9 was first derived by Gross, 
Pisarski and Yaffe and Weiss for SU(N) gauge bosons and for massless fermions p-pHj. 

Each of the four terms deserves comment. The first is the blackbody free energy for two degrees of freedom, and 
depends only on the temperature and the angle 9. The second term, which is the leading correction due to the mass 
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M, often appears in discussions of symmetry restoration at high temperatures with 8 — 0. For example, suppose we 
are calculating the effective potential for a complex scalar field <E>. The mass M 2 is given by the second derivative 
of the classical potential, d 2 U/d§*d$, and depends on the expectation value of the field <I>. If U has the form 



U 



m 2 $*$ + A ($*$) z , then M 2 



-4A$*$. For m > 0, the U(l) symmetry is spontaneously broken at low 



temperature. At high temperature, the M 2 / 12/3 2 term generates a positive mass for the $ field of order T, restoring 
the symmetry. The third term in Vb (6) is linear in T, and non-analytic in M 2 for 9 = 0. It is closely associated 
with the n = Matsubara mode, which is the most infrared singular contribution to a finite temperature functional 
determinant. This term is responsible for non-analytic behavior in finite temperature perturbation theory via the 
summation of ring diagrams. For example, in a scalar theory it gives rise to the A 3 / 2 contribution to the free energy; 
in QED, the contribution is e 3 Q Note how subtractions occur in the I =^ parts of this term to keep these parts 
subleading. The last term is logarithmic in the dimensionless combination [3M and independent of 9. In calculations 
of effective potentials, it typically combines with zero-temperature logarithms in such a way that the temperature T 
sets the scale of running coupling constants at high T. 

From the basic result for Vb, we can build other results. Consider a complex scalar field in the fundamental 
representation of SU (N) . The partition function in a constant background Polyakov loop is given by 
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where we assume all #'s are chosen to lie between and 2tt in accordance with the convention for 9+. 
expression holds for bosons in the adjoint representation. 
For fermions, we have similarly 
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with now used. For fermions in the fundamental representation of SU(N), we may write 
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where the angles 9j must now be chosen to lie between —n and it. 
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V. CONCLUSIONS 



We have found complete, simple expressions for Vg (9) and Vf (9) in the high-temperature limit which generalize 
previously known expressions. Not only are the expressions simple, their derivation is direct and relatively elementary. 
Our formulae reflect in a direct way periodicity in 9, a property which is lost when analytically continuing power 
series in f3/j, to /3fj, + i9. 

As a practical matter, it is natural to ask how accurate both the low- and high-temperature expansions are. The 
low temperature expansion for Vb (#) is an infinite series in n; using the first 10 terms in the series gives an accuracy 
better than 1 part in 10 3 over the entire range to 2ir for temperatures T < 0.25M. The high temperature expansion 
also involves an infinite series, in the the parameter I. In comparison, the high temperature expansion is within 5% 
of the exact answer for all values of 9 at T = 0.5M when terms up to |Z| = 10 are included. The accuracy improves 
substantially as T increases. Both expansions are more accurate when restricted to 9 — 0. 

Our primary interest in these results lies in their application to the study of systems where a non-trivial Polyakov 
loop is expected to occur. The foremost physical example is QCD at finite temperature. The high-temperature form 
of Vb{9) suggests that the Bernoulli polynomials appear naturally in the free energy of SU(N) gauge theories with 
a non-trivial Polyakov loop, essentially as polynomials in the Polyakov loop eigenvalues. In our recent work with 
Miller |jl2), we have used this observation to construct a phenomenological free energy for the quark-gluon plasma 
which reproduces much of the thermodynamic observed in lattice simulations. One can also apply the results obtained 
here to the Savvidy model at finite temperature [JL5| , [L6| . Savvidy originally proposed a model of the QCD vacuum 
in which gluons moved in a constant chromomagnetic field [jl7| . Using low-temperature expansions, we have shown 
that a confining state, where the Polyakov loop expectation value is zero, minimizes the free energy. Using the Bessel 
function identities proven here, we have recently |p_8[] developed a high-temperature expansion for this model which 
shows that a non-zero Polyakov line is favored at high temperature. 



APPENDIX A: DERIVATION OF LOW TEMPERATURE EXPANSIONS 



We begin by expanding the logarithms and performing the angular integrations: 
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This can in turn be reduced using standard recursion relations for modified Bessel functions: 
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APPENDIX B: PROOF OF A BESSEL FUNCTION IDENTITY 



In this appendix, we prove the Bessel function identity 
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We introduce a regulating mass fj,, which will be taken to zero at the end of the calculation, obtaining 
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The first integral gives a sum of screened Coulomb, or Yukawa, potentials 
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and both terms appear to be problematic as /i — > 0. The first term can be made finite in this limit by subtracting the 

contribution at r = <p = for n ^= 0. Using the notation Y] n to denote a summation over all n with the omission of 
the singular term when n — 0, we have 
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The second term in brackets can be evaluated by summing the series 
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and using the identity 
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In the limit \i — > 0, the second term in brackets gives 
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